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Abstract. The SARS epidemic of 2002-3 led to the study of epidemic mod- 
els including management measures and other generalizations of the original 
1927 epidemic model of Kermack and McKendrick. We consider some natural 
extensions of the Kermack-McKendrick model and show that they share the 
main properties of the original model. 



1. Introduction. Almost since the beginning of recorded history there have been 
epidemics. An epidemic may be described as a sudden outbreak of a disease that 
infects a substantial portion of the population in a region before it disappears. In 
the nineteenth century recurrent waves of cholera killed millions in India. The 
influenza epidemic of 1918-1919 killed at least 20 million people overall, more than 
half a million in the United States. One of the questions that first attracted the 
attention of scientists interested in the study of the spread of communicable diseases 
was why diseases would suddenly develop in a community and then disappear just 
as suddenly without infecting the entire community. 

One of the early triumphs of mathematical epidemiology was the formulation of 
a simple model that predicted just such behavior. Kermack and McKendrick [18] 
formulated a model with the population under study being divided into compart- 
ments, namely a susceptible class S, an infective class /, and a removed class R. A 
special case of their model has become known as the Kermack-McKendrick model. 
In this model R is determined once S and / are known, and thus we can drop the 
R equation from the model, leaving the system of two equations: 



The system (1) has been used successfully to fit data from many epidemics, and 
it is very easily analyzed qualitatively. This analysis has been described in many 
references, including [6]. 

The model is based on the following assumptions: 

i. An average member of the population makes contact sufficient to transmit 
infection with /3N others per unit time, where N represents total population 
size. 

ii. A fraction a of infectives leave the infective class per unit time. 
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iii. There is no entry into or departure from the population, except possibly 
through death from the disease. 

Since the fraction of contacts made by an infective that are with a susceptible and 
may therefore produce a new infective is S/N, the rate of new infections is 

PNp = PSI. 

If there are disease deaths, N is a function of t. Let us think of a population of 
size K into which a small number of infectives is introduced, so that S(0) w K, 
1(0) w 0, N(0) = K. 

The model (1) has the following two basic properties shared by more complicated 
epidemic models. 

I. There is a basic reproduction number Rq = /3K/a such that if Ro < 1 there 
is no epidemic, while if Rq > 1 there is an epidemic. 

II. An epidemic eventually dies out, leaving part of the population untouched. 

By an epidemic we mean a situation in which I first increases to a maximum 
and then decreases to zero. The other possibility is that / decreases monotonically 
to zero (no epidemic). The definition of the basic reproduction number Rq is that 
the basic reproduction number is the number of secondary infections caused by a 
single infective introduced into a wholly susceptible population of size K over the 
course of the infection of this single infective. In this situation, an infective makes 
0K contacts in unit time, all with susceptibles and producing new infections, and 
the mean infective period is 1/a; thus, the basic reproduction number is [3K/a. 

To establish the first property, we find orbits of the system (1) in the (S, I) plane 
by integrating 

dl _ I' _ {PS -a)I a_ 

dS~ S> ~ -0SI + P~S' 
Every orbit is a curve in the (S,I) plane, terminating at a point (5oo,0), with 
5oo > and satisfying 

A--|logS = Soo-^IogSoo, (2) 

the final size equation. From this we deduce that lim^oo S(t) = > 0. To show 
that J(t) — > 0, we observe that S + I is a non-negative monotone nonincreasing 
function, which has a limit as t — > oo. Since (S + I) is a smooth function, its 
derivative must approach zero, and this shows that lim^oo I(t) = 0. 

Initially, the number of infectives grows exponentially if and only if Rq > 1, 
because the equation for / may be approximated for t — > by 

I' = (f3K - a)I, 

and the initial exponential growth rate is 

r = (3K-a = a(Ro-l). 

This initial growth rate r may be determined experimentally when an epidemic 
begins. Then, since K and a may be measured, (3 may be calculated as 

However, because of incomplete data and underreporting of cases, this estimate may 
not be very accurate. This inaccuracy is even more pronounced for an outbreak of 
a previously unknown disease, where early cases arc likely to be misdiagnosed. 
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Until the SARS epidemic of 2002-3, most of the work in mathematical epi- 
demiology concentrated on studies of specific diseases or on the interplay between 
epidemiological and demographic effects. There were studies of epidemic models, 
including [4, 5, 6, 13, 14, 15, 16, 22, 24], which concentrated on analysis of the 
course of an epidemic but did not examine the effects of control measures. Studies 
of the SARS epidemic, such as the dynamic models in [2, 9], began with the original 
so-called Kermack-McKendrick model (1) but included control measures to make 
the models more realistic. The purpose of this work is to show that a large vari- 
ety of general epidemic models, including those which have come out of the SARS 
epidemic, have the basic properties I and II. 

2. General contact rates. The assumption in the model (1) of a rate of contacts 
per infective proportional to population size TV, called mass action incidence, was 
used in all the early epidemic models. It is more realistic to assume a contact rate 
that is a nonincreasing function of total population size. For example, a situation in 
which the number of contacts per infective in unit time is constant, called standard 
incidence, is probably a more accurate description for sexually transmitted diseases. 

We generalize the model (1) by replacing assumption (i) with the assumption 
that an average infective makes C(N) contacts in unit time with C'(N) > 0, [1, 7], 
and we define 

«*> - 

It is reasonable to assume f3'(N) < to express the idea of saturation in the number 
of contacts. Then bilinear incidence corresponds to the choice C(N) = (3N, 0(N) — 
(3, and standard incidence corresponds to the choice C(N) = X,/3(N) = A/TV. The 
assumptions C(N) = N(3(N),C'(N) > imply that 

0(N) + Nf3'(N) > 0. (3) 

Some disease transmission models have used a Michaelis-Menten type of inter- 
action of the form 

<™ - 

[7]. Another form based on a mechanistic derivation for pair formation [12] leads 
to an expression of the form 

C(N) = aN / 

l + bN+ VI + 2WV 

Data for diseases transmitted by contact in cities of moderate size [21] suggests a 
good fit is obtained using the form 

C(N) = \N a , 

with a = 0.05. All of these forms satisfy the conditions C'(N) > 0,0' (N) < 0. 
Additional remarks on contact modeling may be found in [19]. If there are no disease 
deaths, so that the total population size is constant, all incidence assumptions 
involving total population size are equivalent, but if there are disease deaths the 
behavior of a model might depend on the form of the incidence. 

An epidemic model in which the incidence is assumed to depend on total pop- 
ulation size must include an equation for total population size. This forces us to 
distinguish between members of the population who die of the disease and members 
of the population who recover with immunity against reinfection. Such a distinc- 
tion was not needed in the original Kermack-McKendrick model, even in its general 
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form. We assume that a fraction / of the members leaving the infective class at 
time t recover and the remaining fraction (1 — /) die of disease. 

To generalize (1) to include a general contact rate, we use S, 7, and N as variables, 
with R= N — S — I. It will be convenient to use N as a model variable in place of R 
because the rate of new infections is now (3(N)SI. We obtain a three-dimensional 
model: 

5" = -0(N)SI 

I' = j3(N)SI-aI (4) 
N' = -(l-f)al. 

We wish to show that model (4) has the same qualitative behavior as model (1), 
namely that it has the properties I and II. 

For model (4), the basic reproduction number is given by 

Kq — 

a 

because a single infective introduced into a wholly susceptible population makes 
C(K) = K(3(K) contacts in unit time, all of which are with susceptibles and thus 
produce new infections, and the mean infective period is 1/a. In addition to the 
basic reproduction number R , there is also a time-dependent running reproduction 
number, that we call R* , representing the number of secondary infections caused 
by a single individual in the population who becomes infective at time t. In this 
situation, an infective makes C(N) — Nf3(N) contacts in unit time, of which a 
fraction S/N are with susceptibles and thus produce new infections. Then it is 
easy to see that for model (4) the running reproduction number is given by 

a 

If R* < 1 for all large t, the epidemic will die out. We may calculate the rate of 
change of the running reproduction number with respect to time, using (3) and (4) 
to find that 

d R * = S'(t)(i(N) + S(t)0'(N)N'(t) = (-P(N)) 2 SI - Sa(l - f)f3'{N) 
dt a a 



< /3(N)SI 



a 



Thus, f t R* < if Nf3(N) > a(l - /), or R* > (1 - f)S/N. This means that R* 
decreases whenever R* > 1. Thus, if 72* < 1 for t = T, then R* < 1 for t > T. 
If i? > 1, then I'(0) — a(R a — 1)7(0) > 0, and an epidemic begins. However, R* 
decreases until it is less than 1 and then remains less than 1. Thus, the epidemic 
will die out. If R < 1, then 7'(0) = a(R - 1)7(0) < 0, 7?* < 1 for all t, and there 
is no epidemic. 

If we use the same approach as was used for (1) to show that S 1 ^ > 0, we obtain 

dl a 

dS = ~ + S(3(N) ' 

and we are unable to proceed because of the dependence on N. However, since 

(3(N) < (3(0), 
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we have 

dI > 1 a 



dS ~ S/?(0) ' 

We may now use the following comparison theorem established by R. Conti [3]. 
Theorem 1. Let y(x), z(x) be the respective solutions of the initial value problems 

y' = f(x,y), y(x )=y 

z ' = g(x,y), z(xo) = yo 

and suppose that g(x,y) < f(x,y) for all x,y. Then 

z{x) < y(x), x > x 
z{x) > y(x), x < xo- 

If we apply this comparison theorem to the solutions I(S) of 

dl „ a 
-1 + 



dS Sf3(N) 
and I (S) of 

dl „ a 

= -1 + 



dS S(3(0) ' 

we see that for S < K we have I(S) < Io(S). Since we have seen in our study of 
(1) that the the graph of Iq(S) does not reach S = 0, it follows that the same is 
true for the graph of I(S), and this shows that > 0. We should note that we 
are assuming here that (3(0) is finite. 

If (3(N) — > oo as N — > 0, we must use a different approach to analyze the limiting 
behavior. If / = 1, the total population size remains equal to the constant K, and 
model (4) reduces to the simpler model (1) with f3 replaced by the constant (3(K). 
If / < 1, then TV is a decreasing function of t, which may be inverted to give t as 
a decreasing function of N, and then we may consider S as a function of N. This 
function satisfies the separable differential equation 

dS_ _ £T_ _ (3(N)S 

dN~ iV 7 ~ a(l - /) ■ ( ' 

We may solve (5) with initial condition S = N = K corresponding to t = by 
separation of variables, obtaining 

t-K 



S(t) = Kexp 



■SZ w (3(N)dN 
a(l-f) 



(6) 



We deduce from (6) that = is possible only if TV — > and 0(N)dN diverges. 
If / > 0, it is clear that it is not possible to have N — > 0, since N(i) > R(t) and 
R(t) takes positive values. Thus = is possible only if / = 0. 

Our conclusion is that for model (4), some susceptibles escape the epidemic, 
as with the simpler model (1), unless all infectives die of disease and /3(N)dN 
diverges. This divergence condition means that (3(N) must be unbounded as N — > 0, 
a condition that is biologically unreasonable. In particular, standard incidence is 
not realistic for small population sizes. A more realistic assumption would be that 
the number of contacts per infective in unit time is linear for small population size 
and saturates for larger population sizes, which rules out the possibility that the 
epidemic could sweep through the entire population [19]. 
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In many infectious diseases, there is an exposed period after the transmission of 
infection from susceptibles to potentially infective members but before these poten- 
tial infectives can transmit infection. A generalization of the epidemic model (4) 
that incorporates an exponentially distributed exposed period with compartments 
S, E, I, R, and total population size N = S + E + I + R is 

S' = -0(N)SI 

E' = (3(N)SI-kE (7) 
I 1 = kE- al 
N' = -(1 -/)«/. 
We obtain, much as in the analysis of (1), 

d(E + I) a 
dS + S0(N)' 

and we may deduce from this that 5^ > 0, as with (4). It follows that (7) and (4) 
have the same asymptotic behavior. 

Some diseases have an asymptomatic stage in which there is some infectivity 
rather than an exposed period. This may be modeled by assuming infectivity 
reduced by a factor e A during an exposed stage. The analogue of model (7) with 
infectivity during the exposed stage and a density-dependent contact rate is 

S' = -l3{N)S{I + e A A) 

A' = (3(N)S(I + e A A) - kA (8) 

/' = kA- al 

N' = -(l-f)al. 

Here, we have used A, for asymptomatic, in place of E as one of the compartments. 
For this model 

K(3{K) K(3{K) 

H a = h £a ■ 

a k 

If we attempt to use the approach used when there is no infectivity during the 
exposed stage, we obtain 

d _ al 

dl { + >~~ /3(N)(I + e A A)S' 

which we are unable to integrate. However, we can use an alternate approach which 
actually comes from the original Kcrmack-McKcndrick analysis. 
We begin by writing 

/' = -S' -A' - al, 1(0) = 
and solving this linear first order differential equation to give 

I(t) = - [\s'( S )+A'( S )}e- a(t - s Us. 
Ja 

A similar integration of the initial value problem 

A' = -kA - S', £7(0) = 

gives 

A{t) = - [ S'{s)e-< t - s) ds. 
Jo 
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Then we have 

s'(t) 

S(t) 

which we integrate with respect to t from to oo, obtaining 



-0(N(t)) [ [S'(s) + A'(s)]e-^ t - s Us + pe A f [-S'(a)]e- K ^ds, 
Jo Jo 



log— = / 0{N{t)) [-S , {s)][e- a <t-') + e A e-« t -*]d8dt 

boo Jo Jo 

-0{N{t)) f f A'(s)e- a ^- s Usdt 
Jo Jo 



10 Jo 



/•OO l>00 

< /3(0) J J [-S\s)][e- a ^ + CAe-^-^dtds 

/'OO pOO 

-j3{K) / / A'{s)e- a{t - s Utds 
Jo Jt 

/•oo />oo 

]/ / [S'{s)][e- a ^ +e A e- KU ]duds 
Jo Jo 



0(0) 



p oo poo 

—0(K) / / A'(s)e- au duds. 
Jo Jo 



Since J °° A'(s)ds = 0, this reduces to 

log^<0(O)(K~S oo )^±^. (9) 
boo na 

This is the analogue of the limit equation, (2), for the simple model, (1), with mass 
action incidence, and shows that Soo > 0. 

3. Models incorporating public-health responses. An actual epidemic differs 
considerably from the idealized models, (1) or (4), as was shown by the SARS 
epidemic of 2002-3. Some notable differences include the following: 

1. Diagnosed infectives may be hospitalized, both for treatment and to isolate 
them from the rest of the population. 

2. Contact tracing of diagnosed infectives may identify people at risk of becoming 
infective, who may be quarantined (instructed to remain at home and avoid 
contacts) and monitored so that they may be isolated immediately if and 
when they become infective. 

3. Isolation may be imperfect; in-hospital transmission of infection was a major 
problem in the SARS epidemic. 

In the SARS epidemic of 2002-3 in-hospital transmission of disease from patients 
to healthcare workers or visitors because of imperfect isolation accounted for many 
of the cases. This points to an essential heterogeneity in disease transmission that 
must be included whenever there is any risk of such transmission. 

All these generalizations have been considered in studies of the SARS epidemic 
of 2002-3. While the ideas were suggested in SARS modelling, they are in fact 
relevant to any epidemic. One beneficial effect of the SARS epidemic has been to 
draw attention to epidemic modeling which may be of great value in coping with 
future epidemics. 

When no vaccine is available, isolation and quarantine are the main measures 
available for attempting to manage an outbreak of a new disease. We assume that 
an epidemic has started, but that the number of infectives is small and almost 
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all members of the population are still susceptible, and we formulate a model to 
describe the course of an epidemic when management measures are begun under 
the following assumptions: 

1 . Asymptomatic members may be infective with infectivity reduced by a factor 
ca,0 < e A < 1. 

2. Asymptomatic members who are not isolated become infective at rate K\. 

3. We introduce a class Q of quarantined members and a class J of isolated 
members. 

4. Asymptomatic members are quarantined at a rate jiA (in practice, a quar- 
antine will also be applied to many susceptiblcs, but we ignore this in the 
model). The effect of this assumption is that some susceptibles make fewer 
contacts than the model assumes. Quarantine is not perfect, but reduces the 
contact rate by a factor eg. 

5. There may be transmission of disease by isolated members, with an infectivity 
factor of ej. 

6. Infectives are diagnosed and isolated at a rate 72/. In addition, quarantined 
members are monitored and isolated immediately when they develop symp- 
toms, at a rate k 2 Q- 

7. Infectives leave the infective class at rate ail, and a fraction /1 of these 
recover, and isolated members leave the isolated class at rate a 2 J , with a 
fraction f 2 recovering. 

These assumptions lead to the SAQIJR model 

S' = -P{N)S[e A A + e A e Q Q + I+ejJ} 
A' = P(N)S[e A A + e A e Q Q + I + ej J] -(« 1+7l )A 
Q' = 1iA-k 2 Q (10) 
I' = k x A - (ai + 72)/ 
J' = K2Q + J2I-U2J 
N' = -(l-/i)a 1 /-(l-/ 2 )a 2 J. 
The model before management measures are begun is (8), the special case 
7i = 72 = «2 = a 2 = h =0, Ki = k, ai=a 

of (10). The model (10) is equivalent to the SARS model of [9] except for the 
extension to a general contact rate in place of standard incidence and the omission 
of immigration and natural death rate terms. The model of [2] is closely related to 
model (10), having isolation but not quarantine, but also distinguishing two classes 
of susceptibles with different susceptibilities. 

We define the control reproduction number R c to be the number of secondary 
infections caused by a single infective in a population consisting essentially only 
of susceptibles with the control measures in place. It is analogous to the basic 
reproduction number, but instead of describing the very beginning of the disease 
outbreak, it describes the beginning of the recognition of the epidemic. We assume 
that this occurs soon enough so that the total population size is still approximately 
K. The basic reproduction number is the value of the control reproduction number 
before management measures are implemented. 

In addition, there is a time-dependent effective reproduction number R* that con- 
tinues to track the number of secondary infections caused by a single infective as the 
epidemic continues with management measures (quarantine of asymptomatics and 
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isolation of symptomatics) in place. It is not difficult to show that if the inflow into 
the population from travellers and new births is small (i.e., if the epidemiological 
time scale is much faster than the demographic time scale) , our model implies that 
R* will become and remain less than unity, so that the epidemic will always die out. 
Even if R c > 1, the epidemic will abate eventually when the effective reproduction 
number becomes less than unity. However, it should be remembered that if the 
epidemic takes so long to die out that there are enough new births and immigrants 
to keep R* > 1, there will be an endemic equilibrium, meaning that the disease will 
establish itself and remain in the population. 

We may calculate R c in the same way as we calculate i? but using the full model 
with quarantined and isolated classes. We obtain 

R ^ = e A K(i(K) | K/3(K) 

R e Q e A Kp(K) | K0(K) Kl | € Q K0{K) lx | ejK0{K)n ll2 | ejK/3(K) 7l 

Di DiD 2 Dik 2 a 2 DiD 2 a 2 D\ 



where D± = 71 + Ki, D 2 =72 + 0:1. 

The running reproduction number R* is the control reproduction number with 
K replaced by N to reflect the change in total population size and multiplied by 
S/N to reflect the fact that the fraction of contacts by an infected member which 
are with a susceptible and thus can produce a new infection is S/N . Thus 



R* = 



e Q e A N/3(N) N0(N) Kl , e Q Nf3{N) ll ejN0(N) 



= Sf3(N)[ 



+ 



D X D 2 + 



DlK 2 

ej«i72 



+ 



a 2 DiD 2 
ej7i 1 



(K172 + JlD 2 ) 



S 
N 



Di DiD 2 D\K 2 a 2 DiD 2 a 2 Di 

Each term of R c has an epidemiological interpretation. The mean duration in A is 
1/D\ with contact rate ca/3, giving a contribution to R c of eEK(3(K)/D\. A frac- 
tion K\jD\ goes from A to /, with contact rate /3 and mean duration l/D 2 , giving a 
contribution of K(3{K)n\/ D\D 2 . A fraction 71/-D1 goes from A to Q, with contact 
rate €q€q(3 and mean duration 1/k 2 , giving a contribution of cqK [3{K)^ij 'Dik 2 . 
A fraction K172/ D\D 2 goes from A to / to J, with a contact rate of ej/3 and a 
mean duration of 1/0:2, giving a contribution of ejK f3(K)nij 2 /a 2 DiD 2 . Finally, 
a fraction 71 /D\ goes from E to Q to J with a contact rate of ejf3 and a mean 
duration of l/a 2 , giving a contribution of ejK (3(K)ji/ 'D\a 2 . The sum of these 
individual contributions gives R c . 

The linearization of (10) at the disease-free equilibrium (K , 0, 0, 0, 0, K) has the 
4x4 matrix 



~e A K(3{K) - ( Kl 

7i 
«i 




■7i) 



e A e Q K(3(K) 


K 2 



K0(K) 


-(01 + 72) 

72 



ejK(3{K) 



— OL 2 



obtained by dropping the zero rows and columns corresponding to S and N. The 
corresponding characteristic equation is a fourth-degree polynomial equation whose 
constant term is positive if R c < 1 and negative if R c > 1. Thus, if R c > 1, there is a 
positive eigenvalue, corresponding to an initial exponential growth rate of solutions 
of (10). If R c < 1, it is possible to show that all eigenvalues of the coefficient matrix 
have negative real part, and thus solutions of (10) die out exponentially [23]. 
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An argument similar to the one used for (4) but technically more complicated 
may be used to show that > for the treatment model (10). If ea = — that 
is, if there is no infectivity in the exposed or asymptomatic stage, — we can give a 
simpler argument, as follows. 

Let U = A + Q + I + J. Then we may consider U as a function of S, and 

dU t a\I + a 2 J 

- 1 



dS /3{N)S{I + ejJ)' 

If we take ci = min(ai,a 2 /ej) we obtain 

dU a 
^ — 1 



dS ~ f3(N)S' 

and we obtain > by the same analysis as was used for (4), provided (3(0) is 
finite. If (3(N) — > oo as N — > 0, we may also use the analysis of (4) since 

dS Sj3(N) 
dN ~ c 2 

if c 2 < mm(ai(l - fi),a 2 (l - fij/tj)- 

Our conclusion is that the asymptotic behavior of the treatment model (10) is 
the same as that of the simpler model (4) . If the control reproduction number R c 
is less than 1, the disease dies out and if R c > 1, there is an epidemic that will die 
out leaving some members of the population untouched. 

To relate and TVqo, we begin by integrating the equations for S + A, Q, I, J, 
and N of (10) with respect to t from t = to t = oo, using the initial conditions 

S(0) + A(0) = N(0) = K, Q(0) = 1(0) = J(0) = 0. 

We obtain 

/•OO 

K-Soo = (ki+7i) / A(s)ds 
Jo 

/•oo />oo 

7i / A(s)ds — k 2 / Q(s)ds 
Jo Jo 

/•oo />oo 

«i / A(s)ds = (ai + 72) / I(s)ds 
Jo Jo 

/•oo /*oo />oo 

k 2 / Q(s)ds = 012 J(s)ds - 72 / I(s)ds 
Jo Jo Jo 

/•OO POO 

K - Noa = (l-h) ai I(s)ds + (1- f 2 )a 2 J(s)ds. 

Jo Jo 

Now we need to express / °° I(s)ds and / °° J(s)ds in terms of / °° A(s)ds. From 
the above relations for integrals we obtain 



Thus we have 
K-N^ 



(ai+72) / I(s)ds — k 1 A(s)ds 
Jo Jo 

f°° t( \ j 7i"i + 7i72 + ^i72 f°° . , 

a 2 / J(s)ds = / A(s)ds. 

Jo "l + 72 Jo 

(1 - /l)«lKl + (1 - /2)(7l a l + 7172 + Kl72) 



Ql + 72 

(1 - /i)ai«i + (1 - /2)(7i Q; i + 7i72 + KiTg) 
(«i +7i)(«i +72) 



/•oo 

[i^-^oo]- (11) 
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This has the form 

K-N oa =c[K- S^} (12) 
with c, the disease death rate, given by 

_ (1 - /QaiKi + (1 ~ /gXTigi + 7i 72 + «i7g) 
(«i + 7i)(«i +72) 

The mean disease death rate may be measured, and this expression gives informa- 
tion about some of the parameters in the model that cannot be measured directly. 
It is easy to see that < c < 1 with c = if and only if fi = f 2 = 1, that is, if 
and only if there are no disease deaths, and c = 1 if and only if /1 = fa = 0, that 
is, if and only if the disease is universally fatal. The argument used to establish 
(11) becomes considerably simpler when applied to model (4) and gives (12) with 
c = 1 - / 

The computations in this section show that the asymptotic behavior of the treat- 
ment model (10) is the same as that of the simpler model (4). If the control repro- 
duction number R c is less than 1, the disease dies out, and if R c > 1, there is an 
epidemic that will die out, leaving some members of the population untouched. 

4. A vaccination model. If a vaccine is available for a disease that threatens an 
epidemic outbreak, a vaccinated class that is protected at least partially against 
infection should be included in a model. While this is probably not relevant for an 
outbreak of a new disease, it would be an important aspect for modeling influenza 
outbreaks or a bioterrorist outbreak of smallpox. 

If there is an outbreak of a disease for which a vaccine is available, then it 
is natural to include the effect of vaccination in an epidemic model. To model 
(4), we add the assumption that in unit time a fraction <f> of the susceptible class is 
vaccinated. The vaccination may reduce but not completely eliminate susceptibility 
to infection. We model this by including a factor a, < o < 1, in the infection 
rate of vaccinated members, with a — meaning that the vaccine is perfectly 
effective and a = 1 meaning that the vaccine has no effect. We assume also that 
the vaccination loses effect at a proportional rate 9. We describe the new model by 
including a vaccinated class V, with 

5" = -p(N)SI -<j>S + 6V 

V = c/>S-a{3(N)VI-8V (13) 

I 1 = f3(N)(S + aV)I-aI 

N' = -(1 - f)al. 

We may view model (13) as a management model for an epidemic modeled ini- 
tially by (4) that is different from model (10). We think of introducing a single 
infective into a population of total size K consisting of 5* susceptibles and V vac- 
cinated members. This infective makes [3(K)(S + aV) contacts in unit time, and 
the mean infective period is 1/a. Thus the corresponding reproduction number is 

L(S(0)+aV(0)). (14) 
a 

There are two essentially different scenarios. The first is an outbreak of a new 
disease for which a vaccine (presumably developed originally for some other dis- 
ease) is available. Then the population would not have been vaccinated before the 
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beginning of the disease outbreak and we would take 5(0) w K, V(0) = 0. Then 
the reproduction number given by (14) would be 

K/3(K) 

Ho = • 

a 

The second scenario is a disease against which the population has been prevac- 
cinated; in this scenario we would assume the population size has reached the 
disease- free equilibrium of (13); namely, 

6 + 4> 6 + <f> 

Then the control reproduction number given by (13) is 

KP(K) 9 + acf, 
a ' 9 + (j> ■ 

From the equation for / in (13), we see that initially / grows exponentially, and 
there is an epidemic if and only if 

0(K)(S(O) + aV(0)) > a, 

that is, if and only if the reproduction number is greater than 1. 

The next step in the analysis of model (13) is to integrate the equations 

(S + V + I)' = -aI 

and 

N' = -(l-f)aI 

to obtain 

K-N 00 = (l-f)[K-S 00 -V 00 ]. 
Next we introduce a new variable T = S+V and then, using S + crV < T, we have 

dl _ a 

-1 + 



dT p{N){S + aV) 

a 

> -1 



P(0)T 

The comparison theorem used previously shows that T is less than the solution 

of 

dl , a 



dT /?(0)T' 
and thus we have > if (3(0) is finite. 

If j3(N) — > oo as N —>■ 0, we may follow the analysis of (7) with T in place of 
S to conclude that the quantity T = S + V , representing the number of members 
of the population not touched by the disease, remains positive unless / = and 
J„ /3(N)dN diverges. Thus the vaccination model (13) retains the basic properties 
(i) and (ii) of the simple epidemic model (1). In particular, in an epidemic model 
with vaccination, it is not possible to have a backward bifurcation. This contrasts 
with the situation in models with births and deaths in which vaccination may lead 
to endemic equilibria when Rq < 1 and backward bifurcations (multiple endemic 
equilibria or endemic equilibria) when R < 1 [8, 10, 11, 20]. 

If a vaccine is available for a particular disease, it would be straightforward 
to formulate a model incorporating both a quarantine and isolation program and 
vaccination. The same approach may be used to establish the basic properties (i) 
and (ii). A natural question, which could be answered using such a model, would be 
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whether vaccination, quarantine and isolation, or a combination of both measures 
would be the most effective management strategy. Presumably, the answer would 
depend on the parameters for the specific disease. 

5. Models incorporating responses by individuals to an epidemic. An 

aspect of epidemics not addressed in any of the models in earlier sections is the 
possibility that some members may reduce the number of contacts they make with 
other members of the population. One way this may occur is if infectives are sick 
enough to be unable to continue their normal activities, even if they do not require 
hospitalization. Another possibility is that if it is known that there is a disease 
outbreak, some susceptibles will voluntarily reduce their contacts and take hygienic 
measures to decrease the danger of becoming infected. In fact, there is reason to 
believe that educational initiatives to encourage safer behavior may be a useful 
public health measure for coping with a disease outbreak [17]. 

We will formulate a model that generalizes (4) and (7) by incorporating both 
a decrease in the number of contacts made by infectives because of a reduced 
level of activity and a decrease in the number of contacts made by susceptibles in 
response to their awareness of infection. We begin with the reduction of contacts 
by infectives. We replace the assumption that all members of the population make 
C(N) contacts in unit time by the assumption that infectives make qC(N) contacts 
in unit time while other members of the population make C(N) contacts in unit 
time, where < q < 1. Then the total number of contacts in unit time is (N — 
I)C(N) + qIC(N), of which a fraction 

N-(l-q)I 

is with an infective. This implies that the number of new infectives in unit time is 

qSIC(N) 
N-(l-q)r 

To model the reduction of contacts by susceptibles we replace the contact rate 
C(N) by a function C(N,I), satisfying the condition 

in addition to the analogues of the conditions imposed earlier on C(N), namely, 

The resulting generalization of model (7) is 

S' = -6{N,I)SI 

E' = 6(N,I)SI-kE (15) 
/' = kE- al 
N' = -(l-f)al 



with 



6(N,Q= qC{N ^ 



N-(l-q)f 

System (15) is identical to the system (7), with (3(N) replaced by 6(N,I), and has 

r> qC(K,o) 

Ho = . 

a 
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To show that Soo > for model (15), we note that 0(N,I) is a nondecreasing 
function of q for all N, I and thus takes its maximum value C(N, I)/N when q = 1. 
Since C{N, I) < C{N, 0) we have 

and now we may apply the comparison theorem as with (7) to show that Soo > 0. 
Thus the qualitative behavior of model (15) is the same as that of (7), and model 
(15) has the basic properties I and II. 

6. Discussion. We have established that general epidemic models behave in the 
same way asymptotically in the sense that there is a basic reproduction number 
which determines whether there will be an epidemic and that an epidemic will 
pass through a population leaving some members untouched. We conjecture that 
this would remain true for more complicated models with more compartments and 
more stages, including models with heterogeneity of mixing. Of course, our models 
assume that the course of the epidemic is rapid enough that demographic effects 
may be ignored. If this is not true, then it would be possible for a disease to become 
endemic. 

The underlying assumptions in all the models we have described are that the sizes 
of the compartments arc large enough that deterministic models arc appropriate 
and that the mixing of members is homogeneous. While these assumptions are 
probably reasonable once an epidemic is well underway, events at the beginning 
of an epidemic may be quite different. To model events with a small number 
of infectives in a population of susceptibles, we should use a branching process 
stochastic model. 

We have been considering an epidemic in a single location, ignoring travel of 
individuals who may be infective between locations. Modern transportation has 
permitted the rapid transfer of infectious diseases over great distances, and an 
aspect of epidemic management that has become important is the screening of 
travellers who may be infective. Epidemic models which include some movement 
into and out of populations are a natural extension of the models considered here. 
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